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Recent spectroscopic studies at terahertz frequencies for a variety of multiferroics endowed with 
both ferroelectric and magnetic orders have revealed the possible emergence of a new collective 
excitation, frequently referred to as electromagnon. It is magnetic origin, but becomes active in 
response to the electric field component of light. Here we give an overview on our recent advance 
in the terahertz time-domain spectroscopy of electromagnons or electric-dipole active magnetic res- 
onances, focused on perovskite manganites — i?Mn03 (i? denotes rare-earth ions). The respective 
electric and magnetic contributions to the observed magnetic resonance are firmly identified by the 
measurements of the light-polarization dependence using a complete set of the crystal orientations. 
We extract general optical features in a variety of the spin ordered phases, including the A-type an- 
tiferromagnetic, coUinear spin ordered, and ferroelectric be and ah spiral spin ordered phases, which 
are realized by tuning the chemical composition of R, temperature, and external magnetic field. In 
addition to the antiferromagnetic resonances of Mn ions driven by the magnetic field component of 
light, we clarify that the electromagnon appears only for light polarized along the a-axis even in 
the coUinear spin ordered phase and grows in intensity with evolution of the spiral spin order, but 
independent of the direction of the spiral spin plane {be or ab) or equivalently the direction of the 
ferroelectric polarization Ps (fs||c or Ps||a). A possible origin of the observed magnetic resonances 
at terahertz frequencies is discussed by comparing the systematic experimental data presented here 
with theoretical considerations based on Heisenberg model. 

PACS numbers: 300.6495 Spectroscopy, terahertz; 160.4760 Optical properties. 



I. INTRODUCTION 



Recently, there is growing interest in the magnctoelec- 
tric (ME) effect of oxide compounds endowed with both 
ferroelectricity and magnetism, termed multiferroics^. 
The ME effect as a generic characteristic of the multi- 
ferroics is manifested by the magnetic field H control of 
the ferroelectric polarization Pg^ or inversely the elec- 
tric field control of the magnetization'^. Such versatile 
cross-correlation phenomena are highly non-trivial, but 
can host the demanding functionality in the future elec- 
tronics; this is the reason why the study of the multi- 
ferroics in a form of bulk as well as thin film is rapidly 
accelerated^'*'^'^. This recent boom is triggered by the 
discovery of the ferroelectricity in orthorhombically dis- 
torted perovskite manganites, /JMnOa, where R repre- 
sents rare-earth ions such as Gd, Tb, Dy, and their solid 
solutions, by Kimura et al.'^-'^'^. Noticeably, the direc- 
tion of Ps can be flopped from the c- to a-axis by an 



application of the external H along the a- or 6-axis^'^'^. 
Within the framework of Landau theory, the non-zero 
component of the ME tensor can produce the linear 
ME effect, which has been known to emerge in non- 
centrosymmetric magnets such as Cr203, in which both 
space-inversion and time-reversal symmetries are simul- 
taneously broken^'^"'^^. However, the gigantic ME ef- 
fect observed in iJMnOa is better understood in terms 
of the phase transition between the multiferroics states; 
this requires the new mechanism of the ferroelectricity 
in these compounds and a variety of experimental ap- 
proaches have been taken to reveal the multiferroic na- 
ture of i?Mn03i''*'5'6^ 

As a typical example of the ferroelectric behavior in 
i?Mn03, we show in Figs. 1(a) and 1(b) the tempera- 
ture dependence of dielectric constant e at 10 kHz and 
Pg along each crystallographic axis of TbMnOa in zero 
H, respectively^'^'*. Below the Neel transition temper- 
ature Tn of 39 K, where the coUinear spins of Mn ions 
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FIG. 1: (Color online) Temperature dependence of (a) dielec- 
tric constant e at 10 kHz and (b) ferroelectric polarization Pa 
of TbMnOa with each crystallographic axis in zero magnetic 
field^ 



sinusoidally order along the 6-axis [see, the schematic il- 
lustration of the collincar spin order in Fig. 2(c)]; e along 
the a-axis tends to slightly increase and shows a kink at 
28 K. Below 28 K, Pg along the c-axis steeply increases 
and reaches the maximum {^-^ 800 fiC/m^) at the low- 
est temperature. Accordingly, e along the c-axis exhibits 
a sharp peak at 28 K, signaling the ferroelectric phase 
transition. In this ferroelectric phase, the presence of the 
be spiral spin structure, where the spins of Mn ions rotate 
within the be plane with the propagation vector (along 
the 6-axis), was identified by neutron scattering experi- 
ments [see, the schematic illustration of the be spiral spin 
order in Fig. 2(e)]^^. 

As an origin of such a magnetically driven ferroelectric- 
ity in multiferroics, Katsura, Nagaosa, and Balatsky have 
proposed the spin-current model^'^, that Pg can be pro- 
duced by the non-coUincar spin order. In this model, the 
ferroelectricity shows up along the direction perpendicu- 
lar to the wave vector and within the spiral spin plane, 
as given by 



Ps (X dj X {S^ X Sj), 



(1) 



where is a unit vector connecting the nearest-neighbor 
spins. Si and Sj, as shown in Fig. 2(e). According to Eq. 
(1), Ps is expected to emerge along the c-axis in i^MnOs 
with the be spiral spin order, being consistent with the 



experimental observation of TbMnOa shown in Fig. 1. 
Furthermore, the relationship between the ferroelectric- 
ity and the spin structure could be directly revealed in 
the spin-polarized neutron scattering by controlling the 
vector chirality defined by {Si x Sj) of TbMnOa by an 
external electric ficld^'*. Therefore, the observed Pg-flop 
from the c- to a-axis by applying the external H can be 
regarded as the flop of the spiral spin plane from be to ab. 
The presence of the ab spiral spin order with Ps along the 
a-axis was recently confirmed for TbMnOs in H along the 
6-axis^^ and Gdo.yTbo.aMnOa in zero H^^. The equiva- 
lent expression to the spin-current model has also been 
obtained phenomenologically^^ and by considering the 
spin-lattice coupling through the Dzyaloshinski-Moriya 
interaction^®. Recently, ME phase diagrams of i?Mn03 
were theoretically investigated on the basis of numeri- 
cal calculations of a microscopic spin model including 
the anisotropy terms and the Dzyaloshinski-Moriya in- 
teraction, which successfully reproduces the experimen- 
tal phase diagrams including flops, emergence, and dis- 
appearance of Ps^^- 

The ME coupling between ferroelectricity and mag- 
netism would also produce the intriguing lower-lying 
spin excitation^°. Since 1960s, the presence of such 
an excitation was theoretically inferred in the ME 
materials^^. This new elementary excitation appears as 
a magnetic resonance in the dielectric constant spectrum 
e(uj) as a response to the electric field component of 
light thus can be termed the electric-dipole active 
magnetic resonance. This is in contrast to the case 
of the magnetic resonance appearing in the magnetic 
permeability spectrum fJ,{Lu) as a result of the spin 
excitation driven by the magnetic field component of 
light i/"^^. Recently, Pimenov et al. have measured 
the optical spectra of TbMnOs and GdMnOs in the 
energy range of 0.4-4.8 meV by using a backward wave 
oscillator (BWO) as a light source, combined with 
Mach-Zehnder interferometry, in which the complex 
optical constants can be extracted in the quasi-optical 
setup^"^. They found a single peak-structure in the 
imaginary part of e{oj), at 2.9 meV for TbMnOs and 
2.5 meV for GdMnOs around 10 K in zero H. These 
absorptions were shown to be allowed when E'^ was 
set parallel to the a-axis. Furthermore, the dramatic 
reduction of the real part of e(uj) as a result of the trans- 
formation from the be spiral spin order to the ^-type 
AFM order, was observed by applying the external H 
along the c-axis. Based on these facts, they claimed the 
possible emergence of the electric-dipole active magnetic 
resonances in these compounds, now frequently refereed 
to as electromagnons. This pioneering work stimulates 
the variety of investigations concerning the low-energy 
spin dynamics of multiferroics such as pcrovskite 

^j^j^Q^23, 24, 25, 26, 27, 28, 29, 30, 31, 32, 33, 34, 35, 36, 37, 38, 39, 40, 41 

hexagonal YMnOg^^ BiFe0343,44^ ^^^^ i?Mn20545.46,47^ 
by using far-infrared optical, Raman scattering, and 
inelastic neutron scattering spectroscopies. 

Within the framework of the spin-current model that 
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can explain the ferroelectric properties observed in 
i^MnOa, Katsura, Balatsky, and N. Nagaosa showed that 
E'^ is possible to drive the oscillation of the spiral spin 
plane in i^MnOs, which produces the magnetic resonance 
at terahertz frequencies along the direction perpendicu- 
lar to the spiral spin plane^^; E'^\\a and E'^\\c in the 
be and ab spiral spin ordered phases, respectively. This 
collective mode can be termed the rotation mode of the 
spiral spin plane, as the example of the manifestation 
of the dynamical ME coupling in multiferroics. Subse- 
quently, Senff et al. performed inelastic neutron scatter- 
ing experiments of TbMnOa in the be spiral spin ordered 
phase and found the magnetic excitation around 2 meV 
at fc = 0^^, the peak position of which is nearly iden- 
tical to that of the electromagnon in TbMnOs revealed 
by BWO spectroscopy^^. Furthermore, they also clari- 
fied that the observed mode is polarized along the a-axis 
perpendicular to the spiral spin plane^^, as predicted by 
theoretical considerations. Therefore, the electromagnon 
observed in i?Mn03 at terahertz frequencies was antici- 
pated to be ascribed to the rotation mode of the spiral 
spin plane. 

Here we show our latest advance in the terahertz 
time-domain spectroscopy on electromagnons or electric- 
dipole active magnetic resonances in multiferroic per- 
ovskite manganites EMnOa^^'^^'^^'^^'^^. Contrary to 
early assignments described above, we provide the com- 
pelling evidence that the electromagnon observed in 
/JMnOa cannot be ascribed to the rotation mode of the 
spiral spin plane; it appears only for i?'^||a, irrespective 
of the direction of the spiral spin plane {he or ab) . This 
conclusion is unambiguously revealed by the systematic 
optical investigations on the spin excitation at terahertz 
frequencies in i?Mn03 based on light-polarization, tem- 
perature, and external H dependence. We clarified the 
full spectral shape of the electromagnon, which consists 
of two peak-structures around 2 mcV and 5-8 meV. The 
lower-lying electromagnon survives even in the coUincar 
spin ordered phase. Below Tn, we also identify the 
presence of the antiferromagnetic resonances of Mn ions 
around 2 meV driven by H'^\\a or H'^\\e, the peak posi- 
tion of which is nearly identical to that of the lower-lying 
electromagnon for £'"||a. 

The format of this article is as follows. In Sec. II, we 
briefly describe the experimental setup for the terahertz 
time-domain spectroscopy we used here and the estimate 
procedure of the complex optical constants from the raw 
data in time domain. Sec. Ill is devoted to show the 
results of the optical spectra at terahertz frequencies for 
a variety of spin ordered phases of i?Mn03, as tuned 
by the ionic radius of i?, temperature, and external H . 
After the brief descriptions of the overall optical spec- 
trum of TbMnOa (Sec. Ill A) and the basic feature of 
i?Mn03 in terms of Mn-O-Mn bond angle (Sec. IIIB), 
we sum up the general optical properties in a variety of 
spin ordered phases by taking DyMnOa (Sec. IIIC)^^ 
and Gdo.rTbo.aMnOa (Sec. IIID)^^ as the examples. 
Wc discuss in Sec. IV the origin of the electromagnons 



in i?Mn03 with theoretical considerations based on the 
Heisenberg model'^^. Summary is given in Sec. V. 



II. METHODS 



We used the standard experimental setup for the ter- 
ahertz time-domain spectroscopy in transmission geom- 
etry. Femtosecond laser pulses delivered from the mode- 
locked Ti:sapphire laser with the center wavelength of 800 
nm, the pulse width of 100 fs, and the repetition rate of 
80 MHz were divided into pump and trigger pulses. The 
pump pulses were irradiated to the ZnTc crystal or to 
the photoswitching device made on the low-tempcrature- 
grown GaAs (LT-GaAs) coupled with the bow-tie an- 
tenna. The detector was the another LT-GaAs coupled 
with the dipole antenna. Although the lower limit of 
the available energy range is restricted by the size of 
each sample, we cover the energy range of 2-9 meV and 
0.8-6 meV by using ZnTe and LT-GaAs terahertz emit- 
ters, respectively. The radiated terahertz pulse was colli- 
mated and focused on the sample by a pair of the off-axis 
paraboloidal mirrors. For the light-polarization depen- 
dence, the wire-grid polarizer was inserted in between 
the off-axis paraboloidal mirrors. 

Single-crystalline samples were grown by the floating- 
zone method*. Specimens with wide ac, ab, and he faces 
were cut from the bowl and each crystallographic axis 
was determined by back Laue photographs. The ob- 
tained specimens were characterized by x-ray diffraction, 
e at 10 kHz, Pg, and magnetization measurements, which 
were all consistent with previous reports in Rcfs.^'^'*. For 
transmission experiments, we polished the specimens to 
the thickness of 100-850 /xm. We carefully confirmed that 
there was no effect of the polishing procedure on the op- 
tical properties at terahertz frequencies of /JMnOa. 

We estimated the optical constants n of /JMnOa with- 
out the Kramers-Kronig transformation. In i?Mn03, 
there are spin excitations driven by both E'^ and H'^ 
in the measured energy range, as we could clarify their 
contributions to fi mainly based on the measurements 
of the complete set of the light-polarization dependence 
(Sec. HI). Due to the emergence of the magnetic reso- 
nances driven by iJ'^, n should be precisely expressed as 
fi = Y^e/i (whereas fi ~ ^/e for the case of non-magnets). 
Wc confirmed that the contribution of ^ is negligible to 
the complex transmission coefficient by the numerical cal- 
culation and thus the effect of H'^ was taken into account 
by adopting h = ^JT\x. Therefore, we used the quantity 
of EyU in this paper. Further details of our estimate pro- 
cedure and the validity of this approach can be found in 

Rcf.29. 
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FIG. 2; (Color online) (a) Magnetoelectric phase diagrams of Tt^MnOs with change of the Mn-O-Mn bond angle, as reproduced 
from Refs.^^'^^ {R represents rare-earth ions). The solid lines are merely the guide to the eyes, (b) Magnified view of (a), 
showing the emergence of the ferroelectric phases for R = Gd, Tb, and Dy. By changing the Mn-O-Mn bond angle, iiMnOs 
shows the large variations of the spin structures with temperature, including the A-type (layer- type) antiferromagnetic (AFM), 
coUinear spin ordered, be and ab spiral spin ordered, and i?-type AFM phases. The schematic illustrations of these phases are 
shown in (c)-(g); Mn ions and their spins are highlighted by circles and arrows, respectively, e^j is the unit vector connecting 
adjacent spins. Si and Sj. According to the spin-current mechanism, as formulated by Eq. (1), the ferroelectric polarization 
Ps appears along the c- and a-axes in the be and ab spiral spin ordered phases, respectively. 



III. EXPERIMENTAL RESULTS 

A. Overall optical spectrum of TbMnOa from 
terahertz to ultraviolet frequencies 

Fh'st, we show the overall optical spectrum of TbMnOa 
in the ferroelectric be spiral spin ordered phase, rang- 
ing from terahertz to ultraviolet frequencies'^^. Figures 
3(a) and 3(b) show the real Re[e/i] and imaginary Im[e/i] 
parts of the e/i spectra of TbMnOs up to 90 meV, re- 
spectively, measured at 12 K. was set parallel to the 
a-axis. Although there is no contribution of fi to the e/i 
spectrum above 8 meV, we used the notation of Re[e/i] 
and Im[e//]. We obtained the spectrum above 10 meV 
by using the Kramers-Kronig transformation; the polar- 
ized reflectivity spectrum was measured in the energy 
ranges of 0.01-0.8 eV and 0.6-36 eV by using Fourier 
transform infrared spectrometer and grating monochro- 
mator, respectively. In the energy range of 10-22 meV, 



we performed both the transmission and reflectance mea- 
surements and thus directly estimated the e/j, spectrum 
without the Kramers-Kronig transformation. Below 10 
meV, terahertz time-domain spectroscopy was used in 
transmission geometry. The inset of Fig. 3(b) shows 
the Im[e/i] spectrum of TbMnOa up to 10 eV, measured 
at 10 K, where the electronic transitions arc dominant. 
The optical transition around 2 cV across the charge- 
transfer gap is clearly identified. The peak-structures 
around 5 eV and 9 eV can be assigned to the transi- 
tions from O 2p to Mn 3d and Tb 5c? states, respectively, 
according to the systematic optical study of the tran- 
sition metal oxides with perovskite structure"*^*'^^ . Al- 
though the GdFcOa-type distortion results in the split- 
ting and mixing of the phonon modes, we can roughly 
classify the character of the observed phonon modes into 
three types — stretching, bending, and external modes, 
which are typical for ideal cubic perovskite structure'*^ 
The optical phonon around 70 meV is assigned to the 
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FIG. 3; (Color online) Overall optical spectrum of TbMnOs, ranging from 2 meV to 10 eV, in the ferroelectric be spiral spin 
ordered phase, measured around 10 K. was set parallel to the a-axis. (a) Real Re[e/i] and (b) imaginary Im[e/i] parts 
of the e/j. spectrum up to 90 meV. The low-energy part of the Im[e/i] spectrum below 21 meV is multiplied by 10. Inset 
shows the Im[e^] spectrum up to 10 eV, measured at 10 K. Below 10 meV, both Re[e/i] and Im[e/x] spectra were obtained by 
terahertz time-domain spectroscopy. Above 10 meV, we used the Fourier transform spectrometer (0.01-0.8 eV) and the grating 
spectrometer (0.3-36 eV). 



stretching mode of Mn-O-Mn. Several modes discerned 
from 25 meV to 60 meV mainly come from the bending 
modes of Mn-O-Mn. The optical phonon modes below 
25 meV can be ascribed to the external modes, which 
correspond to the vibrations of Tb ions. Below 20 meV, 
other contributions to the far-infrared spectrum arc iden- 
tified as four pronounced absorption peak-structures at 
3 meV, 7 meV, 14 meV, and 17 meV, as clearly seen in 
the Im[e^] spectrum, multiphed by 10 [Fig. 3(b)]. The 
14 mcV peak-structure corresponds to the lowest-lying 
external phonon mode, which shows an appreciable cou- 
pling with the lower- lying electromagnon^^'"^-'^'^^. Accord- 
ing to the inelastic neutron scattering experiments^^'^^, 
the upper edge of the magnon band of TbMnOa is lo- 
cated around 8 meV. Therefore, we previously assigned 
the peak-structure at 17 meV to the upper edge of the 
two-magnon band, as the energy of this peak-structure 
just corresponds to the twice of the magnon energy at 
zone edge"^^. However, we recently clarified that this 
absorption peak is diminished in other RMnO^ such as 
DyMn03^^ GdMn03^^ and Eui-^Y^MnOg (x repre- 
sents the nominal composition)^^, although two other 
peak-structures below 10 meV are discerned. There- 
fore, the peak-structure at 17 meV can be ascribed to 
the crystal field excitation of Tb'^"'' ions, contrary to 
our previous assignments^. Further detailed studies con- 
cerning the phonon modes in i?Mn03 will be presented 



elsewhere^^'S^. In the following subsections, we focus on 
optical properties of i?Mn03 below 10 meV. 



B. Overview of terahertz spectra of iiMnOa 

Here we overview the optical spectra at terahertz fre- 
quencies in thermally induced spin ordered phases — A- 
type antiferromagnetic (AFM) and be spiral spin ordered 
phases, realized in a family of -RMnOa (i? = Gd, Tb, 
and Dy). In i?Mn03, the decrease of the ionic radius 
of R or cquivalcntly the decrease of the Mn-O-Mn bond 
angle </>, destabilizes the A- type AFM order^^, in which 
the Mn spins antiferromagnetically stack along the c-axis 
and ferromagnetically order along the o- and 6-axes, as 
schematically shown in Fig. 2(f). The ME phase di- 
agrams of i?Mn03 are reproduced in Figs. 2(a) and 
2(b)^^'^s. LaMnOa (0 = 155.1°) has one eg electron per 
Mn-site (d^) and is an A-type AFM insulator below Tn of 
140 K. With decreasing the ionic radius of R from R = La 
{(j) = 155.1°) to i? = Gd ((/) = 146.2°), Tn dramatically 
decreases from 140 K to 42 K [Fig. 2(a)] as a result of 
the spin frustration caused by competing spin-exchange 
interactions^^. In between i? = Tb (</> = 145.4°) and 
R = Dy {(j) = 144.7°), the ferroelectricity emerges along 
the c-axis, accompanied by the be spiral spin order [Fig. 
2(e)]2.7,8^ as can be clearly seen in the magnified view 
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FIG. 4: (Color online) Low-energy electrodynamics of the spin excitation of (a) GdMnOs at 12 K, (b) TbMnOa at 12 K, and 
(c) DyMnOs at 9 K, for and -ff'^||c. Upper and lower panels show the real Re[e^] and imaginary Im[e/i] parts of the efj, 

spectra, respectively. Around 10 K, the be spiral spin ordered phase is realized for TbMnOa and DyMnOa, while the A-type 
antiferromagnetic phase for GdMnOg. 



[Fig. 2(b)]. At 10 K, GdMnOa and Gdo.rTbo.sMnOs are 
A-typc AFM insulators. Noticeably, around cj) of 146.0°, 
as can be realized in the mixed valence Gdi-^Tb^^MnOa 
in between the paraelectric ^-type AFM GdMnOs and 
the ferroelectric TbMnOa^"^, the ab spiral spin ordered 
phase with Ps||a [Fig. 2(d)] emerges in a narrow range 
of temperature; for example, -Ps||a shows up between 16 
K and 24 K for Gdo.TTbo.gMnOa^^ Below 16 K, the 
A-type AFM phase is stable. For further decreasing the 
ionic radius of R, i?-type AFM order is developed for 
R = Ho [Fig. 2(g)], in which the Mn spins order in the 
sequences of |m along each crystallographic direction. 

We show the Re[e/Lt] and Im[e/i] spectra in upper and 
lower panels of Fig. 4, respectively, for (a) GdMnOa at 
12 K, (b) TbMnOa at 12 K, and (c) DyMnOg at 9 K. 
In these measurements, E'^ and H'^ were set parallel to 
a- and c-axes, respectively, with use of the crystal plates 
with a widest ac face. As can be clearly seen, there is a 
large variation of the optical spectra with the ionic radius 
of R. In the ^-type AFM phase of GdMnOa (0 = 146°) 
at 12 K [Fig. 4(a)], a clear sharp peak-structure can 
be discerned at 2.3 meV in the Im[e^] spectrum, yield- 
ing the maximum magnitude of Im[e//] ~ 1. In accord 
with this, there is a dispersive structure in the Re[e^] 
spectrum. This tiny absorption can be assigned to the 
spin excitation driven by H'^\\c, the details of which are 
discussed in Sec. HID by adopting Gdo.rTbo.sMnOs. 
With decreasing the ionic radius of i?, we can see the 



dramatic modification of the e/i spectra. In the be spi- 
ral spin ordered phase (Ps||c) of TbMnOs at 12 K [Fig. 
4(b)], the magnitude of Im[e/i] increases, which forms 
the pronounced broad continuum-like absorption com- 
posed of two peak-structures at 2.9 meV and 7.4 meV. 
The magnitudes of Im[e/z] for the lower- and higher- lying 
peak-structures reach about 5 and 7, respectively. At 
these peak positions, we can identify the clear disper- 
sive structures in the Re[e/z] spectrum. In DyMnOa with 
the further decreased the ionic radius of R but similar 
with be spiral spin order, the lower-lying peak-structure 
grows in intensity at 9 K and the magnitude of Im[e/.i] 
reaches the maximum about 10 [Fig. 4(c)]. The position 
of the lower-lying peak-structure shifts from 2.9 meV to 
2 meV. On the contrary, the magnitude of Im[e//] of the 
higher-lying peak-structure decreases, but the peak po- 
sition also shifts to the lower-energy from 7.4 meV to 5 
meV; such tendency is also discerned in the Re[e/i] spec- 
trum. The observed remarkable absorptions in TbMnOa 
and DyMnOs can be ascribed to the electromagnons, as 
discussed in detail in Sec. Ill C. Further systematic inves- 
tigations on the variation of the optical spectra in terms 
of R can be found in Ref.^^. In the following subsections, 
we focus on the optical spectra at terahertz frequencies 
in the A-type AFM, be spiral, and ab spiral spin ordered 
phases of i?Mn03, tuned by temperature and external 
H, by taking DyMnOg (Sec. Ill C) and Gdo.yTbo.sMnOg 
(Sec. HID) as examples. 
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FIG. 5: (Color online) Temperature 
dependence of (a) dielectric constant 
e at 10 kHz and (b) ferroelectric po- 
larization Ps of DyMnOs with each 
crystallographic axis in zero magnetic 
field'''^ The magnetic field (|| b) 
dependence of (c) dielectric constant 
along the a-axis at 10 kHz and (d) Pa 
(Psikand PJla). 



C. DyMnOs with be spiral spin order 



In DyMnOa , the collinear spin order of Mn ions evolves 
along the 6-axis below Tn of 39 K. The modulation 
wavevector of Mn ions g^" along the fe-axis is incom- 
mensurate with q^" ~ 0.36''. For decreasing temper- 
ature, the ferroelectricity appears along the c-axis as a 
result of the evolution of the be spiral spin order; the 
temperature dependence of e at 10 kHz and Ps along the 
each crystallographic axis of DyMnOa in zero H are pre- 
sented in Figs. 5(a) and 5(b), respectively^'*. Below the 
ferroelectric transition temperature Tc of 19 K, Pg along 
the c-axis steeply increases up to ~ 2000 ^G/tc? [Fig. 
5(b)]. In accord with this, e along the c-axis exhibits the 
sharp peak at Tc [Fig. 5(a)]. Contrary to the behavior 
of e along the a-axis of TbMnOa with the same he spi- 
ral spin order [Fig. 1(a)], e along the a-axis of DyMnOa 
strongly enhances below Tn and yields the large value 
of e 50) at T^. e along the a-axis is roughly by a 
factor of 2 larger than e along other axes. We show in 
Fig. 5(d) the H dependence of Ps along the a- and c- 
axes, measured at 9 K. By an application of external H 
along the &-axis, Pg along the c-axis is dramatically sup- 
pressed. On the other hand, P, along the a-axis steeply 
increases, accompanied by the huge change of e along the 
a-axis [Fig. 5(c)], when the external H exceeds the crit- 
ical value He of ~ 20 kOe. Namely, the direction of Ps 
can be flopped from the c- to a-axis above or equiva- 
lently the spiral spin plane changes from he to ah. Among 
a family of PMnOa, DyMnOa shows the largest ME ef- 
fect, as exemplified by the remarkable change (300% at 
10 kHz) of e along the a-axis upon the Ps-flop [Fig. 5(c)] 
and by the large Ps [Fig. 5(d)]^'*. Therefore, the tera- 
hertz time-domain spectroscopy of DyMnOa provides the 



useful insights into the basic characteristics of the elec- 
tromagnons observed in PMnOa and their possible role 
in the ME effect. 

First, we clarify the selection-rule of the spin excita- 
tion in the be spiral spin ordered phase of DyMnOs based 
on the light-polarization dependence (for both i?" and 
i/") using the complete set of the crystal surface plates 
(ac, ah, and 6c)'^^. These measurements are indispensable 
to experimentally distinguish the respective electric and 
magnetic contributions to the e/i spectrum. Figures 6(a) 
and 6(b) present the Re[e/i] and Im[e^] spectra in the be 
spiral spin ordered phase, respectively, measured around 
10 K. There is a remarkable optical anisotropy with re- 
spect to i?" and H'^ . According to the inelastic neu- 
tron scattering experiments for the case of TbMnOa^^'^^, 
the crystalline-electric-field excitation of / electrons lies 
around 4.7 mcV. However, the observed remarkable opti- 
cal anisotropy as well as the negligible absorption at 11 K 
for ii^'^llc and H'^\\h [open circles in Figs. 6(a) and 6(b)] 
can exclude the possible emergence of such an excitation 
in the c^l spectrum in the measured energy range. 

In the Im[e/i] spectrum for £"^11 a and ff"||c at 9 K 
[closed squares in Figs. 6(a) and 6(b)], there is a pro- 
nounced broad absorption up to 10 meV, which con- 
sists of two peak-structures around 2 mcV and 6 meV. 
Accordingly, the clear dispersive structures are visible 
in the Rc[e^] spectrum. The presence of the lower- 
lying absorption has been reported for TbMnOs at 2.9 
meV23, GdMnOg at 2.5 meV^^, and Eui_:i,Y:j,Mn03 at 
3 meV^^'^^. Among them, the magnitude of Im[e^] of 
DyMnOs reaches about 10, which is roughly by a fac- 
tor of 2-5 larger than that of TbMnOs and GdMnOs, as 
also seen in Fig. 4. The observed remarkable character- 
istics, including the positions for peak-structures and the 
spectral shape, can be also identified for E'^\\a and H'^\\b 
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FIG. 6: (Color online) Light-polarization dependence of the spin excitations in the be spiral spin ordered phase of DyMnOa, 
measured around 10 K, using a complete set of the crystal faces (ac, ab, and be). The crystal orientations with respect to _E" 
and are indicated in the figures. Upper and lower panels show the real Re[e/^] and imaginary Im[e^] parts of the e/x spectra 
(symbols), respectively. Im[e/i] spectra shown in (d) are vertically offset for clarify. Note that the scales of the vertical axes 
in (a) and (b) are different in (c) and (d), respectively. The solid lines shown in (a) and (b) are results of a least-square fit to 
reproduce lower- and higher-lying peak-structures by assuming two Lorentz oscillators for e. On the other hand, the e/x spectra 
shown in (c) and (d) can be reproduced by two Lorentz oscillators for e and p. 



at 10 K [closed circles in Figs. 6(a) and 6(b)]. There- 
fore, the observed broad absorption can be ascribed to 
the electric-dipole active mode only along the a-axis. We 
can reproduce the e/i spectra for E'^\\a by adopting two 
Lorentz oscillators of e for peak-structures, as indicated 
by the solid lines. It yields Re[e^(tt' 0)] of ^ 32, which 
is less than e along the a-axis at 10 kHz (~ 50) [Fig. 5(a)]. 
A noticeable characteristic is the large optical anisotropy 
of the Re[e/i] spectrum. Re[e/i] for iJ'^||a at low-energy 
yields the large value of Re[e/i] of 25-30, which is about 
2 times larger than that along the other axes, i.e., 15.5- 
17.5. This large anisotropy ratio (^ 2) is comparable to 
that of e at 10 kHz (~ 2) [Fig. 5(a)], indicating that the 
anisotropic dielectric response is extended from kilohertz 
to terahertz frequencies associated the gigantic electric- 
dipole active excitation. 

On the contrary to the E'^\\a case, we found tiny ab- 
sorptions around 2.2 meV for -ff"|ja and i?"||c, whose 
peak energy is nearly identical to that of the lower-lying 
peak-structure for i?'^||a. Figures 6(c) and 6(d) show the 
Re[e^] and Im[e/i] spectra, respectively, measured around 
10 K. There is a single peak-structure around 2.2 meV 
in the Im[e^] spectrum [Fig. 6(d)] with a clear disper- 
sive structure in the Re[e/i] spectrum [Fig. 6(c)] for £'"||6 
and -ff^jja at 8 K (closed triangles), although the magni- 
tudes are an order of magnitude smaller than the £'"||a 



peaks. This peak-structure can be assigned to the spin 
excitation driven by -ff'^lja as the nearly identical spec- 
tra signature can be observed for E'^\\c and H'^\\a at 
6 K (open squares), while the absolute value of the e/i 
is slightly different. We also discern the broad peak- 
structure around 2.3 meV for E'^Wb and i?"||c at 9 K 
(crosses). For the case of E'^\\a and -ff"||c, the possible 
peak feature driven by is completely masked by the 
intense electric-dipole active absorption for £^'^||a. This 
peak-structure for E'^\\h and H'^\\c diminishes above Tn, 
as in the same manner to the spin excitation by H^\\a. 
Therefore, this peak-structure is ascribed to the spin ex- 
citation driven by H'^\\c^ details of which are discussed 
in Sec. HID by adopting Gdo.TTbo.sMnOs. 

In the following, we mainly focus on the observed gi- 
gantic absorption for E'^Wa. Figures 7(a) and 7(b) show 
the temperature dependence of the Rc[e/i] and Im[e/i] 
spectra, respectively, for E'^\\a and H'^\\c. At 245 K, 
there is no remarkable absorption. The slight accumula- 
tion of the Im[e/i] spectrum above 6 meV is due to the 
contribution of the optical phonon absorption of the per- 
ovskite structure, located around 14 meV^'* (see, also the 
overall optical spectrum of TbMnOs shown in Fig. 3). 
With decreasing temperature, the magnitude of Im[e/i] 
below 6 meV tends to increase. Around Tn of 42 K, 
Im[e/i] forms the peak-structure around 2 meV. Accord- 
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FIG. 7: (Color online) Temperature dependence of (a) real 
Re[e/i] and (b) imaginary Im[£/i] parts of the selected e/i spec- 
tra of DyMnOa for E'^\\a and -ff^ljc. The solid lines are re- 
sults of a least-square fit to reproduce lower- and higher-lying 
peak-structures in the be spiral spin ordered phase below 19 
K by assuming two Lorentz oscillators for e. The selected efi 
spectra of DyMnOs in the magnetic field applied along the 
&-axis, measured at 7 K, are also included. Above 19 K, the 
data can be fitted by the single Lorentz oscillator. 



ingly, there is the dispersive structure in the Re[e/L<] spec- 
trum. In the coUinear spin ordered phase below Tn, the 
absorption grows in intensity, while the peak position 
shifts to the lower energy. In the be spiral spin ordered 
phase below Tc of 19 K, the absorption becomes promi- 
nent and finally reaches the maximum of Im[e^] 10 
at low temperature. In this ferroelectric phase, we also 
identify the additional peak-structure around 6 meV in 
the Im[e/x] spectrum, which forms the broad continuum- 
like absorption up to 10 meV. Its peak position shifts to 
the lower energy with decreasing temperature, while the 
lower-lying peak position shifts to the higher energy. 

To further clarify the origin of the broad continuum- 
like absorption for E'^\\a, we show in Fig. 8(a) the tem- 
perature dependence of the integrated spectral weight per 
Mn-site (Ncs) defined as, 



iVn- 



2mQV 



wim [e (ti> ) /i ( w ) ] , 



(2) 



where mo is the free electron mass, e the elementary 
charge, and V the unit-cell volume. We chose Wi = 0.8 
meV and W2 = 9.5 mcV to fully cover the broad absorp- 
tion. By the definition in Eq. (2), the estimated A^off in- 
cludes the contribution of /_i for H'^\\c, which is crudely es- 
timated to be as small as 0.02 x 10^^. Above Tn, A'off for 
-E'^lja (circles) shows the negligible temperature depen- 
dence. However, A^cft gradually increases below Tn and 
sharply enhances below Tc [see, also in the inset of Fig. 
8(a)]. Therefore, we can conclude that the observed ab- 
sorption for £^"^110 is a magnetic in origin, namely, i5'^||a 
electromagnon. On the contrary, TVcff for H'^\\a (squares) 
exhibits the slight enhancement below Tn and shows the 
maximum at Tc [Fig. 8(a)]. We also plot the contribu- 
tion of the background absorption (a dotted line) as N^s 
includes the contribution of e for E'^Wc [the inset of Fig. 
8(a)]. 

In i?Mn03, the rotation mode of the spiral spin plane 
was proposed at the early stage as the origin of the elec- 
tromagnon observed at terahertz frequencies^^. Within 
this picture, the rotation mode of the spiral plane has 
the particular selection-rule for E^; it would become ac- 
tive along the a- and c-axes in the be and ab spiral spin 
ordered phases, respectively. Such an electromagnon sce- 
nario was considered to explain the inelastic neutron scat- 
tering spectra of TbMnOs in the be spiral spin ordered 
phase^^. To test the origin of the electromagnon for 
iJ^lla, we studied the effect of on the optical properties 
in DyMnOa for E'^Wa and i?"||c. In DyMnOa, applying 
H along the &-axis can induce the ab spiral spin ordered 
phase with Ps||a; the direction of Pg can be flopped from 
the c- to a-axis at 20 kOe, as presented in Fig. 5(d). 
During the Pg-flop, the incommensurate g^" keeps nearly 
constant (~ 0.37) for DyMnOs^"', contrary to the case of 
TbMnOa, which shows the incommensurate (g^" = 0.38) 
to commensurate (g^" = 1/4) transition^^. 

In Figs. 7(a) and 7(b), we also include the H depen- 
dence of Re[e^] and Im[e^] spectra, respectively, mea- 
sured at 7 K. The external H was applied along the b- 
axis. We confirmed that the Pg-Aop occurs at 20 kOe 
[Fig. 5(d)]. By application of H, the position of the 
lower-lying peak-structure shifts to the lower energy and 
its intensity grows. Noticeably, we can still see the elec- 
tromagnon for i?'^||a showing the gigantic absorption, 
even in the ab spiral spin ordered phase, as exemplified 
by the e/i spectra at 33 kOe and 59 kOe. We plot in 
Fig. 8(b) the N^s (circles) as a function of H up to 70 
kOe. The integrated range of A^cff is chosen as oji ~ 0.7 
mcV and 1^2 = 5 meV. As can be seen, there is negligi- 
ble H effect on N^s even when the ab spiral spin order 
is developed above He of 20 kOe apart from the slight 
enhancement of Ncs- These observations are different 
from the reported result of TbMnOs, in which the lower- 
lying peak-structure disappears by an application of the 
external H^^. In the experiments in Ref.^^, the external 
H was applied to the c-axis to destruct the ferroelectric 
order, accompanied by the transformation from the be 
spiral spin order to the A-type AFM order. According 
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FIG. 8: (Color online) Integrated spectral weight per Mn-site A'^off, as defined by Eq. (2) in the text, of the ac surface crystal 
plate of DyMnOs as a function of (a) temperature in zero magnetic field and (b) magnetic field at 7 K. A^cfi for -E'^Uti and 
H'^Wa are represented by circles and squares, respectively. The solid lines are merely the guide to the eyes. Inset of (a) shows 
the magnified view below 45 K. For comparison, A'cff for -ff"||a are multiplied by 40. The horizontal dashed line in the inset of 
(a), also multiplied by 40, represents the estimated contribution of the background absorption for -E"^||c to TVeff for H'^Wa. The 
vertical solid lines in (a) and the inset of (a) indicate the ferroelectric transition temperature of 19 K and the Neel transition 
temperature of 42 K. In (b), the magnetic field of 20 kOe, at which the direction of the ferroelectric polarization is fiopped 
from the c- to a-axis or the spiral spin plane changes from be to ah, is indicated by the vertical solid line. 



to the ME phase diagrams of TbMnOa and DyMnOa^, 
the application of H along the b- or a-axis is needed to 
induce the Pg-Aop from the c- to a-axis or cquivalently 
the rotation of the spiral spin plane from he to a6, as we 
performed here. To confirm the survival of the electro- 
magnon for -E'^jja, we also measured the efj. spectra for 
E'^\\c and H'^Wa in H along the 6-axis. H dependence 
of A^eff for E'^\\c and H'^\\a (squares) is presented in Fig. 
8(b). We chose the integrated range as Wi = 2.4 meV 
and W2 = 4 meV. Above H^, the electromagnon arising 
from the rotational motion of the ah spiral spin plane, is 
expected to appear for E'^\\c in the ah spiral spin ordered 
phase. However, A^eff for E'^\\c and i/"||a is considerably 
smaller than A^eff for E'^Wa and i?'^||c by a factor of 100, 
and there is no remarkable change of A^eff at H^- 



D. Gdo.rTbo.aMnOa with ab spiral spin order 

Here we present the results of Gdo.yTbo.sMnOa in zero 
H^^ , which represent the generic optical properties of 
the ah spiral spin ordered phase of i?Mn03. Among a 
family of -RMnOa, Gdo.yTbo.aMnOa is a rare example 
as it exhibits the ferroclcctricity along the a-axis even 
in zero H [Fig. 2(b)]. This is in contrast to the cases 
of GdMnOs, TbMnOs, and DyMnOa, where Ps||a phase 
was only induced when the external H is applied along 
the b- or a-axis, as described in Sec. HI C for the case 
of DyMnOs. We show in Fig. 9(a) the temperature de- 
pendences of Ps and e at 10 kHz along the a-axis for 
Gdo.yTbo.sMnOs'^^. With decreasing temperature, the 
paraclectric coUinear spin order evolves below Tn of 42 
K along the 6-axis, as in the similar manner to other 



i^MnOa showing the he spiral spin order in zero H . Below 
Tc of 24 K, Ps steeply increases and e exhibits the sharp 
peak. Finally, Ps reaches the maximum (~ 280 /^C/m^) 
at 17 K. The estimated Ps is comparable to that along 
the c-axis of TbMnOs in zero H [Fig. l(b)]2'^ Below 
16 K, Ps suddenly vanishes upon the development of the 
A-type AFM order. In this A-type AFM phase, the fi- 
nite magnetization emerges along the c-axis as a result of 
the slight canting of the Mn spins. In Gdo.yTbo.sMnOs, 
the presence of the ah spiral spin ordered phase in zero 
H was recently confirmed by the polarized neutron scat- 
tering experiments with use of an isomer of Gd ion to 
prevent the large cross-section of neutron of Gd ions^®. 
Moreover, g^" along the 6-axis was found to be 0.25, 
which is identical to that in the ah spiral spin ordered 
phase of TbMnOs induced by H^^ . In accord with above 
facts, the ah spiral spin ordered phase of TbMnOs in H 
is smoothly connected with that of Gdo.TTbo.aMnOa in 
zero H , as can be seen in the comparison of ME phase 
diagrams for TbMnOa and Gdo.TTbo.sMnOs, shown in 
Figs. 10(a) and 10(b), respectively. Therefore, we can 
study the generic spin excitation at terahertz frequencies 
in the ah spiral spin ordered phase, which provides the 
further insights into the nature of the observed electro- 
magnon for ii^'^lla in i?Mn03. 

First, we show the light-polarization dependence of the 
e^L spectrum for Gdo.yTbo.sMnOs using a complete set of 
the crystal faces to clarify the light-polarization selection- 
rule of the spin excitation in the ah spiral spin ordered 
phase. Figures 11(a) and 11(b) present the Re[e^] and 
Im[e/i] spectra, respectively, measured around 17 K. For 
these measurements, we used the ac, ah, and he faces to 
distinguish the electric and magnetic contributions. In 
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FIG. 9: (Color online) (a) Thermally-induced ferroelectricity 
in the ab spiral spin ordered phase of Gdo.yTbo.sMnOs^^. The 
ferroelectric polarization Pa emerges along the a-axis between 
16 K and 24 K, accompanied by the change of the dielectric 
constant e along the a-axis at 10 kHz. The critical temper- 
atures for ^-type antiferromagnetic (AFM) order with the 
weak ferromagnetism along the c-axis, ab spiral spin order, 
and coUinear spin order, are 16 K, 24 K, and 42 K, respec- 
tively, as indicated by vertical lines, (b) Temperature depen- 
dence of the integrated spectral weight per Mn-site A'cff, as 
defined by Eq. (2) in the text, of the ac surface crystal plate 
of Gdo.yTbo.sMnOs. A'eft for -E"||a and H"||a are represented 
by circles and squares, respectively. For comparison, A'cff for 
H'^\\a are multiplied by 10. The solid lines in (b) are merely 
the guide to the eyes. 



analogy to the spectra of DyMnOa with the be spiral 
spin order (Fig. 6), we observe the broad continuum-like 
absorption at 17 K (closed squares), which consists of two 
peak-structures around 2 meV and 8 meV, when E'^ and 
were set parallel to a- and c-axcs, respectively. Such 
a broad absorption was assigned to the electromagnon in 
the be spiral spin ordered phase of DyMnOa (Sec. Ill C). 
The observed absorption of Gdo.yTbo.sMnOs was also 
assigned to electromagnon because of the closely similar 
spectral feature for E'^ \ \ a and H'^ \ \ b (closed circles) . Even 
though we checked the reproducibly of the spectra by 
changing the thickness of the sample, the slight discrep- 
ancy of the magnitude of Im[e/x] spectrum is identified. 
This is perhaps due to the slight variation of Gd-to-Tb 
ratio as Gdo.yTbo.aMnOs locates the critical region in 
between the paraelectric A-type AFM and ferroelectric 
be spiral spin ordered phases^'^, as presented in the ME 
phase diagram [Fig. 2(b)]. The magnitude of Re[e/i] for 
E'^\\a at 17 K reaches about 24 at 1 meV (Fig. 11), which 
is comparable to e (^ 28) along the a-axis at 10 kHz [Fig. 
9(a)]. The magnitude of Re[e/x] along the a-axis is larger 
than Re[e/x] for other axes by a factor of 1.5, similar to 
the case for DyMnOs with be spiral spin order (Fig. 6). 



A broad single peak-structure is also discerned for 
E'^\\e and iJ'^Ija at 17 K (open squares), where the peak 
energy nearly matches the peak position of the lower- 
lying electromagnon for E'^\\a [Fig. 11(b)]. However, the 
magnitude of lm[e/i] at the peak energy was estimated to 
be 0.4, which is an order of magnitude smaller than that 
for £'"11 a. This absorption can be ascribed to the spin 
excitation driven by H'^\\a because of the similar spectral 
signatures for £''^||6 and fl^'^lja at 21 K (closed triangles). 
Furthermore, we observe the tiny absorption for £''^||6 
and -ff"||c at 20 K (crosses), which becomes prominent 
in the e/x spectrum at 11 K in the A-type AFM phase, as 
presented in Fig. 12(c). Upon the diminishment of the a6 
spiral spin order, responsible for the gigantic contribution 
of e to e/i, we clearly see the single-peak-structure in the 
efi spectrum at 11 K for E'^\\a and -ff"||c. By comparing 
the spectrum at 11 K for £"^||a and H'^\\e with that at 
10 K for £'^116 and H'^We [Fig. 12(c)], this is attributed 
to the spin excitation driven by H'^\\e. 

These spin excitations for 7J'^||a and 7J'^||c can be in- 
terpreted as antiferromagnetic resonances (AFMRs) of 
Mn ions, which are also observed in canted AFM phases 
of Euo.gYo.iMnOg for H'^\\a'^^ and Lai-^rSr^rMnOg {x < 
0.1) for H'^\\c^^. In fact, the fc = magnon in the A-type 
AFM phase of LaMnOa locates around 3 meV^''. There- 
fore, the measured efx spectrum is considered to consist 
of a sharp resonance for n and broad background absorp- 
tion for e, as given by 

Im[e/i] = Re[e]Im[/i] + Im[e]Re[/x], (3) 

which can be phcnomenologically expressed with two 
Lorentz oscillators for e and /i. The solid lines are the 
results of the least-square fit to the data using Eq. (3), 
which can reproduce the efi spectra for Gdo.yTbo.aMnOa 
[Figs. 11(c), 11(d), 12(b), and 12(c)] and DyMnOa [Figs. 
6(c) and 6(d)]. 

To further see the light-polarized spectral change in 
a variety of the spin ordered phases, we plot in upper 
and lower panels of Fig. 12 the temperature variation of 
Rc[e/x] and Im[e/x] spectra, respectively, for (a) i?'^||a and 
i/'^ljc, (b) E'^We and i?'^||a, and (c) £"||6 and iJ"||c, in 
the A-type AFM phase below 16 K, ferroelectric ab spiral 
spin ordered phase with Ps\\a, between 16 K and 24 K, 
paraelectric collinear spin ordered phase between 24 K 
and 42 K, and paraelectric paramagnetic phase above 42 
K. 

At 10 K in the A-type AFM phase for £;'^||c and iJ'^||a 
[Fig. 12(b)], the clear peak-structure is identified at 
2.1 meV, which can be ascribed to the AFMR driven 
by iJ"||a, as describe above. With increasing temper- 
ature, Im[e/i] tends to decrease in magnitude and the 
peak width becomes broad. In addition, the peak posi- 
tion shifts to the lower energy from 2.1 meV at 10 K to 
1.6 meV at 33 K. Above Tn of 42 K, the peak-structure 
is diminished, as can be seen in the efi spectrum at 50 K. 
For the case of E'^Wb and -ff"||c polarization [Fig. 12(c)], 
a similar tendency is observed. 

On the other hand, there is a remarkable temperature 
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FIG. 10: (Color online) Magnetoelectric phase diagrams of 
(a) TbMnOa and (b) Gdo.yTbo.sMnOs. Data indicated by 
closed circles are taken from Refs.^'*'^^ and the solid lines with 
shaded area are merely the guide to the eyes. The be spiral 
spin plane of TbMnOs, acting as a source of the ferroelectric 
polarization Ps along the c-axis, can be flopped to the ab spi- 
ral spin plane (Ps||a) by an application of the magnetic fleld 
H along the ti-ajcis. This magnetically-induced ab spiral spin 
ordered phase smoothly connects with the thermally-induced 
ab spiral spin ordered phase of Gdo.rTbo.sMnOs. The tera- 
hertz measurements in the paraelectric collinear spin ordered 
phase, ab spiral spin ordered phase, and A-type antiferro- 
magnetic (AFM) phase of Gdo.yTbo.sMnOs were performed 
in zero H, as indicated by the right-side vertical arrow in (b). 



variation of the optical spectra for E'^\\a and H'^Wc [Fig. 
12(a)]. In the ^-type AFM phase at 11 K, the sharp 
peak-structure can be discerned at 2.1 meV, ascribed 
to the conventional AFMR by H'^Wc, described above. 
With increasing temperature, the magnitude of Im[e/x] 
dramatically enhances and the electromagnon emerges 
as the broad continuum-like absorption, as exemplified 
by the efj. spectrum at 17 K in the ab spiral spin ordered 
phase. The intensity of the electromagnon absorption 
decreases when the ab spiral spin ordered is diminished 
above 24 K. However, the electromagnon survives even 
in the colhnear spin ordered phase as a broad absorption 
band [lower panel of Fig. 12(a)], the characteristic en- 
ergy of which was estimated to be about 2.9 meV. Above 
Tn of 42 K, the broadened absorption subsists, while the 
intensity of the electromagnon is completely damped, as 
seen in the efj. spectrum at 68 K. 

To be more quantitative, we show in Fig. 9(b) the 
temperature dependence of A'off (squares) for E'^\\a; the 
integrated range for iVcft [see, Eq. (2)] was chosen from 
u!i = 1.8 meV to u!2 = 6.2 meV to characterize the spin 
excitation spectral range. On cooling, iVofi increases be- 
low Tn and sharply enhances below Tc. This behavior 
is in contrast to the behavior of N^s due to the n com- 
ponent AFMR with E'^\\c and H'^\\a (circles); the latter 



is nearly temperature independent above Tn and grad- 
ually increases below Tn- N^s for iJ'^||a dramatically 
suppresses when the ab spiral spin order transforms to 
the A-type AFM order below 16 K, which contributes to 
the reduction of e at 10 kHz [Fig. 9(a)]. 

In early studies^^, the electromagnon for E'^\\a was 
ascribed to the rotation mode of the spiral spin plane. 
In this picture, the electromagnon would become active 
in the ab spiral spin ordered phase of Gdo.yTbo.aMnOa 
when E'^ was set parallel to the c-axis, perpendicular to 
the spiral spin plane. Figures 11(a) and 11(b) present 
Re[e/i] and Im[e/i] spectra (open circles) for £"^||c and 
H'^\\b in the ab spiral spin ordered phase, measured at 16 
K. As can be seen, there is negligible absorption in the 
measured energy range. We also measured the efj. spectra 
down to 0.8 meV with use of the LT-GaAs terahertz emit- 
ters coupled with the bow-tie antenna as a light source. 
However, the same tendency is discerned. The present 
measurements using a complete set of crystal faces can 
clearly exclude the possibility that the observed electro- 
magnon can be ascribed to the rotation mode of the spiral 
spin plane. Combined with the results presented in Sec. 
HI C for DyMnOs in the 7J-induced ab spiral spin state 
(Ps||a), we can firmly conclude that there is a unique 
selection-rule along the a-axis for the electromagnons in 
.RMnOa, irrespective of the direction of the spiral spin 
plane {be or ab). 



IV. DISCUSSION 

Based on the systematic experimental investigations 
with varying the ionic radius of i?, light-polarization, 
temperature, and external H (Sec. HI), we can ex- 
tract the general features of i?Mn03 in a variety of the 
spin ordered phases, as summarized below. There are 
spin excitations driven by E'^ and H'^ at terahertz fre- 
quencies in i?Mn03. In the ferroelectric spiral spin or- 
dered phase, the electromagnon appears only along the 
a-axis. It spreads in the energy range of 1-10 meV as 
a pronounced continuum-like absorption with two peak- 
structures around 2 meV and 5-8 meV, rather than the 
single peak-structure (1-5 meV) as previously reported 
for TbMnOg^^. Accordingly, Re[e/u(w 0)] for E'^ along 
the a-axis is a factor of 2 larger than those for other 
polarization. This large optical anisotropy ~ 2 is com- 
parable to the anisotropy of e at 10 kHz. The lower- 
lying peak-structure of the electromagnon can be iden- 
tified, though broadened, even in the collinear spin or- 
dered phase but suddenly disappears in the A-type AFM 
phase. In addition, we observed weak but sharp single 
peak-structure for -ff^Ha and i?'^||c in the energy range 
of 1-4 meV, whose peak positions are nearly identical to 
that of the lower- lying electromagnon for £'"||a. These 
peak-structures can be assigned to the AFMRs of Mn- 
ions, which are pronounced in the A-type AFM phase. 

In the experimental data presented in Sec. Ill, we 
could reveal the new features of the optical properties of 
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FIG. 11: (Color online) Light-polarization dependence of the spin excitations in the ah spiral spin ordered phase of 
Gdo.yTbo.sMnOs, measured around 17 K, using a complete set of the crystal faces (ac, ah, and 6c). The crystal orienta- 
tions with respect to E'^ and H"^ are indicated in the figures. Upper and lower panels show the real Re[e/i] and imaginary 
Im[e^] parts of the e/i spectra (symbols), respectively. Im[e/i] spectra shown in (d) are vertically offset for clarify. Note that 
the scales of the vertical axes in (a) and (b) are different in (c) and (d), respectively. The solid lines shown in (a) and (b) are 
results of a least-square fit to reproduce lower- and higher-lying peak-structures by assuming two Lorentz oscillators for e. On 
the other hand, the e/x spectra shown in (c) and (d) can be reproduced by two Lorentz oscillators for e and \i. 



iJMnOa. As an origin of the clcctromagnon for E^^a, 
the rotation mode of the spiral spin plane has been 
proposed^^ based on the spin-current mechanism that 
can explain the emergence of ferroelectricity and the Pg- 
flop by the external H in i?Mn03^'^. Actually, this mech- 
anism was considered to explain the inelastic neutron 
scattering spectrum of TbMnOa in the he spiral spin 
ordered phase^^; the one of the low-energy branch of 
the magnon band around 2 meV was assigned to the 
electromagnon. Noticeably, the clcctromagnon model 
predicts the unique selection-rule in terms of the light- 
polarization; it becomes active for E'^^a and E'^^c in the 
he and ah spiral spin ordered phases, respectively. How- 
ever, we have clearly revealed that the electromagnon 
shows up only along the a-axis, irrespective of the di- 
rection of the spiral spin plane (6c or a6), by measur- 
ing the (.[i spectra in the ah spiral spin ordered phases 
in DyMnOs induced by the external H as well as in 
Gdo.rTbo.sMnOa induced by temperature. There is no 
remarkable peak-structure in the e// spectrum for £''^||c 
in the ah spiral spin ordered phase down to 0.8 meV. 
Based on systematic experimental data described above, 
we conclude that the electromagnon in the energy range 
of 1-10 meV cannot be ascribed to the rotation mode of 
the spiral spin plane. Recently, the negligible effect of 
the external H on the electromagnon is also confirmed in 



TbMnOa by applying the external H along the &-axis'^ ■'^^ 
or the a-axis^^, which is consistent with the results of 
DyMnOa in i/^g 

Electromagnon can only contribute to the change of e 
along the a-axis [Fig. 5(a)], as evidenced by the consis- 
tency of the anisotropy ratio of Re[e/x] at terahertz and 
kilohertz frequencies (~ 2). However, there is a still dis- 
crepancy of the absolute values of Re[e/i] at terahertz and 
kilohertz frequencies for the case of DyMnOa, the reason 
of which can be understood by the presence of the ferro- 
electric domain wall that produces the additional contri- 
bution to e, as recently revealed by the dielectric spec- 
troscopy up to 10 MHz of DyMnOs^^; the e spectrum 
shows the relaxation type dispersion with a relaxation 
rate of ^ 10^^ Hz and the magnitude of e at 10 MHz was 
estimated to be ~ 35, which seems to smoothly connect 
with the low-energy part of the measured Re[e/x] spec- 
trum at terahertz frequencies (Re[e/x([j — ^ 0)] ~ 25-30). 

As an origin of the electromagnon for E^\\a, we pro- 
posed the conventional electric-dipole active two-magnon 
excitation due to the exchange-striction mechanism^^, 
since the observed electromagnon is independent of the 
direction of the spiral spin plane and spreads in the wide 
energy range of 1-10 meV despite the fact that the fc = 
magnon driven by -ff" locates around 2 meV^^'^^''^^'^'^. 
Contrary to this scenario, it was found more recently 



14 



Gclo,7Tbo.3Mn03 

(a) £"||a,H"||c 



(b) ElcHla 



,(C) £l/),H"||c 




2 4 6 8 10 

Photon Energy (meV) 



2 4 6 8 10 

Photon Energy (meV) 



2 4 6 8 10 

Photon Energy (meV) 



FIG. 12: (Color online) Temperature dependence of the real Re[e^] (upper panels) and imaginary Im[e^] (lower panels) parts 
of the selected efi spectra of Gdo.rTbo.aMnOa for (a) £"||a and H'^\\c, (b) E'^Wc and H'^\\a, and (c) and H'^\\c. lm[efi] 

spectra shown in (b) and (c) are vertically offset for clarify. Note that the scale of the vertical axis in the upper panel of (a) is 
different in that of (b) and (c). The solid lines shown in (a) are results of a least-square fit to reproduce lower- and higher- lying 
peak-structures by assuming two Lorentz oscillators for e. The e/i spectra shown in (b) and (c) can be reproduced by two 
Lorontz oscillators for e and /i. 



TABLE I: The nearest-neighbor ferromagnetic interaction Ji, 
next-nearest-neighbor antiferromagnetic interaction J2, and 
the interlayer antiferromagnetic interaction Jc for DyMnOs, 
TbMnOa, and Gdo.yTbo.sMnOs^'^'^'''*^ 

Compounds Ji (meV) J2/IJ1I Jc/\Ji\ 



DyMnOa 


-0.71 


1.22 


1.5 


TbMnOa 


-0.81 


0.78 


2.0 


Gdo.rTbo.aMnOa 


-0.82 


0.75 


2.1 



that the symmetric exchange mechanism can generate 
"one-magnon" excitation when the ordered spins are 
non-coUinear, as in the present case of cycloidal spin 
ordcr'^"''^^ . In the following, wc introduce theoretical con- 
siderations on i^MnOa based on Heiscnberg model and 
compare with experimental data presented in Sec. III. 

The magnetic behaviors of the cycloidal ordered mag- 
nets such as i?Mn03 can be described by a three- 
dimensional S = 2 frustrated Heiscnberg model: 



(4) 



r 



i.i. 



where Ji (< 0) is the nearest- neighbor ferromagnetic 
interaction, and J2 and Jc (> 0), are next-nearest- 
neighbor and interlayer antiferromagnetic interactions, 
respectively (see Fig. 13). Due to the frustration between 
Ji and J2, spiral spin states are the ground state when 
the condition J2/\Ji\ > 0.5 is satisfied, where the spiral 
angle 6 is given by cos 6* = — J1/2J2. The respective in- 
teractions for DyMnOs, TbMnOa, and Gdo.yTbo.aMnOa 
were estimated and summarized in Table. l35'37_ jj 



term is a uniaxial anisotropy term, which fixes the di- 
rection of the spiral spin plane, i.e., a = a for DyMnOa 
and TbMnOs to realize a &c-cycloidal state and a = c 
for Gdp.yTbo.sMnOa to stabilize an a6-cycloidal state. 
We assume that = 0.2|Ji|, which reproduces the 
anisotropy gap observed in the magnon dispersions of 
TbMnOa revealed by inelastic neutron scattering exper- 
iments. ^^'^^ 

The lower-lying magnetic excitations of the model on 
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— - 

J3 

Spin 




the basis of Eq. (4) can be described by the linear spin 
wave theory'^^'^". After rotations of local spin axis, we 
make Holstcin-Primakoff approximations as usual, 



'2S 



i=aj, S'f-S'-aJa;. (5) 



Then, the spin Hamiltonian [Eq. (4)] is diagonalized as 



Hq ~ hujqa^jaq + const. 



(6) 



FIG. 13: (Color online) Heisenberg model for RMnOs with 
nearest-neighbor ferromagnetic interactions Ji (< 0), antifer- 
romagnetic interactions for next-nearest-neighbor along the b- 
axis J2 (> 0), and interlayer interactions Jc (> 0). The unit 
cell is written by dot-dashed lines, where four inequivalent 
spin sites exist due to the orthorhombic lattice distortions. 
Thin (blue) arrows describe the ground state spin configura- 
tion for the be cycloidal state. 



In Eq. (6), the spin wave frequencies hujq are given by 



(7) 



where 



Aa 



B„ = 



2Ji cos a — cos - cos cos — — 

2 2 a b 



J2 cos 29 



J: 



'cos^ 



Jl . 2^ T^Qb 

— — sm - cos cos — — 

2 2 a b 



J: 



cos ■ 



2nqh 



Jr COS 



) cos + ,/, + -, 
D 

T' 



(8) 
(9) 



and a spin wave creation aj and annihilation aq operators 
are defined as 

Uq = CqUq + Sqa}_q, (10) 
aLg = CgoLg + Sqaq, (11) 

where Cq and Sq are coefficients, as given by Cq = 
2Bq/{-[2Aq-2{Al-BlY/^Y+^BlY/\^ndsq^{2Aq- 
2{Al - Blf/^)/{-[2Aq - 2[A\ - Blf'^f + 4S2}i/2. 

The calculated magnon dispersions along the 6-axis 
are shown in Fig. 14. Note that the uniaxial anisotropy 
term hardly affects dispersion relations except for q ^ 
(0,g^",l) [g^" EEE e/{TT/b)]. The magnon dispersions 
calculated with the parameters in Table I'^^ can repro- 
duce those observed by inelastic neutron scattering ex- 
periments in TbMnOa^^'^^. The peak positions of mag- 
netic excitations in spiral ordered phase are also esti- 
mated from hujq by considering spin modes for AFMRs. 
The results are summarized in Table II. Here, AFMRs for 
iJ^II cycloidal spin plane (i?" _L cycloidal spin plane) is 
regarded as an excitation of the magnon at g = (0, g^", 0) 



[q — (0, 0, 1)]. Note that the condition H _L cycloidal spin 
plane corresponds to H'^\\a for DyMnOa and TbMnOs, 
and H'^Wc for GdcrTbcaMnOa, and cycloidal spin 
plane to and i?"||c for DyMnOs and TbMnOs, 

and ff"||a and for Gdo.yTbo.sMnOs. These mag- 

netic resonances are observable due to anisotropy like 
Dzyaloshinski-Moriya interactions. The peak positions 
of the AFMRs at terahertz frequencies driven by iJ" 
29,31,32 g^j^g consistent with that estimated from the spin 
wave theory'^^ . In this way, the simple Heisenberg model 
is a good approximation to describe the magnetic proper- 
ties of DyMnOs, TbMnOa, and Gdo.rTbo.sMnOs at low 
temperatures. Note that, to reproduce the whole com- 
plex ME phase diagram of fJMnOa®, further inclusion of 
anisotropy terms, and Dzyaloshinski-Moriya interactions, 
IS necessary . 

Concerning the electromagnons in i?Mn03, it was ex- 
perimentally revealed that at least two modes, lower- 
lying mode around 2 meV and the higher-lying mode 
aroimd 5-8 meV, are induced for E'^\\a (Table II). The 
higher-lying mode can be ascribed to the the zone bound- 
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FIG. 14: (Color online) Magnon dispersions along the 6-axis 
for DyMnOs, TbMnOs, and Gdo.yTbo.sMnOs. See Table I 
for the parameter values, qc represents the wavenumber along 
the c-axis. 



ary mode, i.e., magnon at g = (0,1,0)^'*'"^^. Such a 
magnon is electrically induced through a symmetric term 
of spin dependent polarizations = 11^(5'^ • Sj), which 
vanishes when a center of inversion is located at the mid- 
dle of the bond connecting site i and j. In i?Mn03, 
polarizations are realized due to the 2>x^ — /'i'lp — 
orbital ordering at Mn sites'^^ and/or orthorhombic lat- 
tice distortions^'', as schematically shown in Fig. 15(a). 
Orbital ordering induces the effective polarization along 
a bond, i.e., if = 7r/4, and distortions that perpendic- 
ular to the bond, i.e., ip = —it /A, where ip is the angle 
from the a-axis. Thus, in general, (/3 is a free parame- 
ter in — 7r/4 < ip < 7r/4, however, the selection-rule for 
the one-magnon excitation is independent of ip as we will 
show in the following. Such an effective polarization is 
also obtained by a microscopic theory for the Mn-O-Mn 
bond^"'^. Pbnra symmetry restricts the polarization pat- 
tern, as schematically shown in Fig. 15(a). In such a po- 
larization structure, we can easily see that one-magnon 
can be induced by for spins with the cycloidal spin 
structure only for E'^Wa condition. On the other hand, 
the origin of the lower-lying mode is still an open ques- 
tion. 

In the external electric field E, Hamiltonian Ti can be 
described as 7i = Tig — E- P^. Therefore, a spin structure 
is modified to increase (decrease) the expectation values 
{Si ■ S'j), when E ■ P^ on the bond is positive (negative). 
For the spin structure in RhlnO^, the propagation vec- 
tor of the cycloidal state is aligned to the &-axis and spins 
along the a-axis are uniform. For the case of i5'^||a stim- 
ulation, the modulation of P^ is uniform along the a-axis, 
while that is staggered along the &-axis. Thus, all spins 
are modulated simultaneously, as schematically shown in 
Fig. 15(b), which can produce the effective coupling of 



the spin structures with E'^ . Thus, spins are oscillated 
by the effective transverse staggered field [Fig. 15(b)] by 
E'^\\a, which can induce the electric-dipole active one- 
magnon resonance. On the others hand, for the case 
of £''^11 6 stimulation, the modulation of Pe along the a- 
axis is staggered, where the effective fields cancels out, 
as schematically shown in Fig. 15(c). In this way, mod- 
ulated polarization does not couple with the spin struc- 
tures and, thus no one-magnon resonance occurs. Note 
that these features are independent of the direction of the 
cycloidal spin plane, since P^. (x Si ■ Sj. In this way, the 
observed selection-rule in RMnOs is understood straight- 
forwardly. 

These processes can be realized by representing the 
symmetric spin dependent polarization Pg = Si - Sj by 
using spin wave operators: 

^Pe^nfe6)(«L-a,.J, (12) 

where fl{qzb) ~ {iS\/ SNU cos ip sin 9{cq^^ — Sq^^), 0, 0) 
and qzb = (0,1,0). Using Eq. (12), the imaginary part 
of complex electric polarizability tensor Im^aa at zero 
temperature, which represents an absorption, is obtained 
from Kubo formula: 

Imxaa(w) = NS^U^ cos^ ipsm"^ 9 

X (Cq,6 - Sq^J'^6{uj - UJzb)Sa. a- (13) 

IniXaa('i') has a peak at LUzb due to an electromagnon ab- 
sorption. Such a peak-structure is consistent with the ob- 
served higher- lying mode around 5-8 meV in DyMnOa^^, 
TbMn033\ and Gdo.TTbo.sMnOs^^ on the assumption 
that Im[e/i] ~ Im^aa'^^- 

V. SUMMARY 

On the basis of the measurements of the ionic radius 
of R, light-polarization, temperature, and magnetic field 
dependence, we uncovered the unique optical features of 
the spin excitations at terahertz frequencies in multifer- 
roic perovskite manganites, i^MnOs. We clearly identi- 
fied that the electromagnon appears for E'^\\a, irrespec- 
tive of the direction of the spiral spin plane {be or ah) 
or equivalently irrespective of the direction of the ferro- 
electric polarization (Ps||c or Ps||a); the direct proof is 
provided by the spectroscopic studies on the ab spiral 
spin ordered phases of DyMnOa induced by H along the 
fo-axis as well as Gdo.rTbo.aMnOs induced by tempera- 
ture. The observed electromagnon is broadly distributed 
over the measured energy range of 1-10 meV, which con- 
sists of two peak-structures around 2 meV and 5-8 meV. 
We also identified the fc = AFMRs of Mn ions for i?" || a 
and -ff^llc below Tn in the narrow energy range of 1-4 
meV, which become prominent when the A-type AFM or- 
der evolves. The AFMR of Mn spins appears as the sharp 
peak-structure around 2 meV, which is nearly identical 
to the peak position of the lower-lying peak-structure of 




FIG. 15: (Color online) (a) Cycloidal spin structures, orbital ordering pattern on Mn sites, distortions of O sites, and direction 
of symmetric spin-dependent polarizations Pe- The directions of Pe, which are defined by an angle (p, are shown by thick 
(brown) arrows. Thin (blue) arrows describe the ground state spin configuration for the ab cycloidal state, (b) For the case of 
-E"|ja polarization, the magnitudes of Pe are modified, as shown on each bond due to spin structure modifications. Thin bold 
arrows are spin structures modulated by the electric fields and dotted arrows are spin structures without electric fields. Thus 
oscillating electric field can induce the effective transverse staggered fields shown by (red) arrows, which can act as a source of 
the electric-dipole active one-magnon resonance at zone boundary, (c) For the case of i5'^||b polarization, the phase of i5" ■ Pe 
is staggered, whereas spins along the a-axis are uniform. In this case, the effective fields cancels out, and no one-magnon 
resonance occurs. 



TABLE II: Comparison of the observed peak positions (in a unit of meV) for electromagnons and antiferromagnetic res- 
onances (AFMRs) with theoretical estimations based on Heisenberg model for DyMnOa, TbMnOa, and Gdo.yTbo.sMnOs. 
Electromagnon emerges as the broad continuum-like band with two peak-structures, which become active only along the a-axis 
(i5"||a), while AFMR appears as the single sharp peak-structure for -ff"||a and -ff"||c. 



Experiment Theory 



Resonance 


Electromagnon 


AFMR 


Electromagnon 


AFMR 


Condition 






E-\\a 




DyMnOs 


2.2 5.0 


2.6 2.9 


5.1 


2.4 3.4 


TbMnOa 


2.9 7.4 


3.0 


7.5 


1.8 2.6 


Gdo.yTbo.sMnOa 


2.8 8.0 


2.0 2.0 


8.0 


2.2 1.6 



the electromagnon for E'^\\a. The electromagnon sur- 
vives even in the collinear spin ordered phase, though 
much broadened, above the ferroelectric transition tem- 
perature, but disappears in the ^-type AFM phase. We 
introduce here one of the possible scenarios to explain 
the observed unique light-polarization selection-rule of 
the electromagnon based on the Heisenberg model on the 
spiral spins in the orbital ordered state. With this model, 
the higher-lying electromagnon around 5-8 meV is as- 
signed to the electric-dipole active one-magnon excitation 
at zone boundary, whose peak position agrees with the 
observation. However, there is a still mystery about the 
origin of the lower-lying electromagnon, which becomes 
prominent with decreasing the ionic radius of R. There- 



fore, further theoretical considerations are needed to fully 
explain whole spectral shape of the electromagnons at 
terahertz frequencies. 

Acknowledgments 

We thank J. Fujioka and Y. Ikebe for their support to 
the measurements. We also thank H. Katsura, F. Ka- 
gawa, M. Mochizuki, Y. Taguchi, and N. Nagaosa for 
valuable discussion from an early stage of our study. This 
work was in part supported by Grant-In- Aids for Scien- 
tific Research (16076205 and 20340086) from the Min- 
istry of Education, Culture, Sports, Science and Tech- 
nology (MEXT), Japan. 



* Electronic address: n-kida@erato-mf.t.u-tokyo. ac.jp 



^ Y. Tokura, "Multiferroics as quantum electromagnets," 



18 



Science 312, 1481-1482 (2006). 

T. Kimura, T. Goto, H. Shintani, K. Ishizaka, T. Arima, 
and Y. Tokura, "Magnetic control of ferroelectric polariza- 
tion," Nature (London) 426, 55-58 (2003). 
T. Lottermoser, T. Lonkai, U. Amann, D. Hohlwein, J. 
Ihringer, and M. Fiebig, "Magnetic phase control by an 
electric field," Nature (London) 430, 541-544 (2004). 
S-W. Cheong and M. Mostovoy, "Multiferroics: a magnetic 
twist for ferroelectricity," Nat. Mater. 6, 13-20 (2007). 
R. Ramesh and N. A. Spaldin, "Multiferroics: progress and 
prospects in thin films," Nat. Mater. 6, 21-29 (2007). 
Y. Tokura, "Multiferroics-toward strong coupling between 
magnetism and polarization in a solid," J. Mag. Mag. 
Mater. 310, 1145-1150 (2007). 

T. Goto, T. Kimura, G. Lawes, A. P. Ramirez, and Y. 
Tokura, "Ferroelectricity and giant magnetocapacitance in 
perovskite rare-earth manganites," Phys. Rev. Lett. 92, 
257201 (2004). 

T. Kimura, G. Lawes, T. Goto, Y. Tokura, and A. P. 
Ramirez, "Magnetoelectric phase diagrams of orthorhom- 
bic TiMnOs," Phys. Rev. B 71, 224425 (2005). 
W. Eerenstein, N. D. Mathur, and J. F. Scott, "Multi- 
ferroic and magnetoelectric materials," Nature (London) 
442, 759-765 (2006). 

T. H. O'Dell, The Electrodynamics of Magneto-electric 
Media (North-HoUand, Amsterdam, 1970). 
Magnetoelectric Interaction Phenomena in Crystals (eds. 
A. J. Freeman and H. Schmid) (Gordon and Breach, Lon- 
don, 1975). 

M. Kenzelmann, A. B. Harris, S. Jonas, C. Broholm, J. 
Schefer, S. B. Kim, C. L. Zhang, S.-W. Cheong, O. P. Vajk, 
and J. W. Lynn, "Magnetic inversion symmetry break- 
ing and ferroelectricity in TbMnOa," Phys. Rev. Lett. 95, 
087206 (2005). 

H. Katsura, N. Nagaosa, and A. V. Balatsky, "Spin current 
and magnetoelectric effect in noncollinear magnets," Phys. 
Rev. Lett. 95, 057205 (2005). 

Y. Yamasaki, H. Sagayama, T. Goto, M. Matsuura, K. 
Hirota, T. Arima, and Y. Tokura, "Electric control of spin 
helicity in a magnetic ferroelectric," Phys. Rev. Lett. 98, 
147204 (2007). 

N. Aliouane, K. Schmalzl, D. Senff, A. Maljuk, K. Prokes, 
M. Braden, and D. N. Argyriou, "Magnetic field in- 
duced flop of cycloidal spin order in muliferroic TbMnOs: 
The magnetic structure of the P\\a phase," Preprint at 
http: //arxiv.org/cond-mat/0902.1890 (2009). 
Y. Yamasaki, H. Sagayama, N. Abe, T. Arima, K. Sasai, 
M. Matsuura, K. Hirota, D. Okuyama, Y. Noda, and Y. 
Tokura, "Cycloidal spin order in the a-axis polarized fer- 
roelectric phase of orthorhombic perovskite manganite," 
Phys. Rev. Lett. 101, 097204 (2008). 

M. Mostovoy, "Ferroelectricity in spiral magnets," Phys. 
Rev. Lett. 96, 067601 (2006). 

L A. Sergienko and E. Datotto, "Role of the 
Dzyaloshinskii-Moriya interaction in multiferroic per- 
ovskites," Phys. Rev. B 73, 094434 (2006). 
M. Mochizuki and N. Furukawa, "Mechanism of lattice- 
distortion-induced electric-polarization flop in the mul- 
tiferroic perovskite manganites," accepted for publica- 
tion in J. Phys. Soc. Jpn. May Issue. Preprint at 
http://arXiv:0811.4015 (2008). 

G. A. Smolenski and I. E. Chupis, "Ferroelectromagnets," 
Usp. Fiz. Nauk 137, 415-448 (1982) [Sov. Phys. Usp. 25, 
475-493 (1982)]. 



V. G. Bar'yakhtar and I. E. Chupis, "Quantum theory of 
oscillations in a ferroelectric ferromagnet. Fiz. Tverd. Tela 
(Leningrad) 11, 3242-3247 (1969) [Sov. Phys. Sohd State 
11, 2628-2631 (1970)]. 

A. M. Balbashov, G. V. Kozlov, A. A. Mukhin, and A. 
S. Prokhorov, "Submillimeter Spectroscopy of Antiferro- 
magnetic Dielectrics. Rare-earth Orthoferrites," in High 
Frequency Processes in Magnetic Materials, (eds. G. Srini- 
vasan and A. N. Slavin) pp56-98 (World Scientific, Singa- 
pore, 1995). 

A. Pimenov, A. A. Mukhin, V. Yu. Ivanov, V. D. Travkin, 
A. M. Balbashov, and A. Loidl, "Possible evidence for elec- 
tromagnon in multiferroic manganites," Nat. Phys. 2, 97- 
100 (2006). 

A. Pimenov, T. Rudolf, F. Mayr, A. Loidl, A. A. Mukhin, 
and A. M. Balbashov, "Coupling of phonons and elec- 
tromagnons in GdMnOs," Phys. Rev. B 74, 100403(R) 
(2006). 

H. Katsura, A. V. Balatsky, and N. Nagaosa, "Dynamical 
magnetoelectric coupling in helical magnet," Phys. Rev. 
Lett. 98, 027203 (2007). 

D. Senff, P. Link, K. Hradil, A. Hiess, L. P. Regnault, Y. 
Sidis, N. Aliouane, D. N. Argyriou, and M. Braden, "Mag- 
netic excitations in TbMnOs: Evidence for a hybridized 
soft mode," Phys. Rev. Lett. 98, 137206 (2007). 
R. Valdes Auilar, A. B. Sushkov, C. L. Zhang, Y. J. Choi, 
S.-W. Cheong, and H. D. Drew, "Colossal magnon-phonon 
coupling in multiferroic Euo.75Yo.25Mn03," Phys. Rev. B 
76, 060404(R) (2007). 

A. Pimenov, A. Loidl, A. A. Mukhin, V. D. Travkin, V. 
Yu. Ivanov, and A. M. Balbashov, "Terahertz spectroscopy 
of electromagnons in Eui-ajYi^MnOs," Phys. Rev. B 77, 
014438 (2008). 

N. Kida, Y. Ikebe, Y. Takahashi, J. P. He, Y. Kaneko, 
Y. Yamasaki, R. Shimano, T. Arima, N. Nagaosa, and 
Y. Tokura, "Electrically driven spin excitation in the fer- 
roelectric magnet DyMnOs," Phys. Rev. B 78, 104414 
(2008). 

D. Senff, P. Link, N. Aliouane, D. N. Argyriou, and 
M. Braden, "Field dependence of magnetic correlations 
through the polarization flop transition in multiferroic 
TbMnOa: Evidence for a magnetic memory effect," Phys. 
Rev. B 77, 174419 (2008). 

Y. Takahashi, N. Kida, Y. Yamasaki, J. Fujioka, T. Arima, 
R. Shimano, S. Miyahara, M. Mochizuki, N. Furukawa, and 
Y. Tokura, "Evidence for an electric-dipole active contin- 
uum band of spin excitations in multiferroic TbMnOs," 
Phys. Rev. Lett. 101, 187201 (2008). 

N. Kida, Y. Yamasaki, R. Shimano, T. Arima, and 
Y. Tokura, "Electric-dipole active two-magnon excita- 
tion in ab spiral spin phase of a ferroelectric magnet 
Gdo.TTbo.sMnOs," J. Phys. Soc. Jpn. 77, 123704 (2008). 
N. Kida, Y. Yamasaki, J. P. He, Y. Kaneko, Y. Ikebe, Y. 
Takahashi, R. Shimano, T. Arima, N. Nagaosa, and Y. 
Tokura, "Electrically driven spin excitation at THz fre- 
quencies in be and ab spiral spin phases of perovskite man- 
ganites," J. Phys. Conf. Ser. 148, 012038 (2009); Proc. 
LXIII Yamada Conf. Photo-induced Phase Transition and 
Cooperative Phenomena. 

R. Valdes Auilar, M. Mostovoy, A. B. Sushkov, C. L. 
Zhang, Y. J. Choi, S.-W. Cheong, and H. D. Drew, "Ori- 
gin of electromagnon excitations in RMnOs," Phys. Rev. 
Lett. 102, 047203 (2009). 

S. Miyahara and N. Furukawa, "Theory of electric field 



19 



induced one-magnon resonance in cycloidal spin magnets," 
Preprint at http://arxiv.org/cond-mat/0811.4082 (2008). 
A. Pimenov, A. Shuvaev, A. Loidl, F. Schrettle, A. A. 
Mukhin, V. D. Travkin, V. Yu. Ivanov, and A. M. Bal- 
bashov, " Magnetic and magnetoelectric excitations in 
TbMnOa," Phys. Rev. Lett. 102, 107203 (2009). 
J. S. Lee, N. Kida, S. Miyahara, Y. Yamasaki, Y. Taka- 
hashi, R. Shimano, N. Furukawa, and Y. Tokura, "System- 
atics of electromagnons in the spiral spin-ordered states of 
TiMnOs," accepted for publication in Rapid Communica- 
tion Section of Phys. Rev. B (2009). 

Y. Takahashi, Y. Yamasaki, N. Kida, Y. Kaneko, T. 
Arima, R. Shimano, and Y. Tokura, "Comprehensive study 
on electromagnons and their coupling with optical phonons 
in Eui_:,Y^Mn03 (x = 0.1,0.2,0.3,0.4, and 0.45)," sub- 
mitted for publication (2009). 

J. S. Lee, N. Kida, Y. Yamasaki, M. Mochizuki, R. Shi- 
mano, and Y. Tokura, "Lattice dynamics of multiferroic 
perovskite manganites", manuscript in preparation. 
A. Pimenov, A. M. Shuvaev, A. A. Mukhin, and A. Loidl, 
"Electromagnons in multiferroic manganites," J. Phys. 
Condens. Matter 20, 434209 (2008). 

D. Senff, N. Aliouane, D. N. Argyriou, A. Hiess, L. P. Reg- 
nault, P. Link, K. Hradil, Y. Sidis, and M. Braden, "Mag- 
netic excitations in a cycloidal magnet: the magnon spec- 
trum of multiferroic TbMnOg," J. Phys. Condens. Matter 
20, 434212 (2008). 

S. Pailhes, X. Fabreges, L. P. Regnault, L. Pinsard-Godart, 
I. Mirebeau, F. Moussa, M. Hennion, and S. Petit, "Hy- 
brid Goldstone modes in multiferroic YMnOs studied by 
polarized inelastic neutron scattering," Phys. Rev. B 79, 
134409 (2009). 

M. Cazayous, Y. Gallais, A. Sacuto, R. de Sousa, D. Lebeu- 
gle, and D. Colson, "Possible observation of cycloidal elec- 
tromagnons in BiFeOs," Phys. Rev. Lett. 101, 037601 
(2008). 

M. K. Singh, R. S. Katiyar, and J. F. Scott, "New magnetic 
phase transitions in BiFeOs," J- Phys. Condens. Matter 
20, 252203 (2008). 

E. Golovenchits and V. Sanina, "Magnetic and magneto- 
electric dynamics in RMn205 (R = Gd and Eu)," J. Phys.: 
Condens. Matter 16, 4325-4334 (2004). 

A. B. Sushkov, R. Valdes Auilar, S. Park, S-W. Cheong, 
and H. D. Drew, "Electromagnons in multiferroic YMn205 
and TbMn205," Phys. Rev. Lett. 98, 027202 (2007). 
C. Fang and J. Hu, "An effective model of magnetoelec- 
tricity in multiferroic RMn2 05," EPL 82, 57005 (2008). 
T. Arima, Y. Tokura and J. B. Torrance, "Variation of 
optical gap in perovskite-type 3d transition-metal oxides," 
Phys. Rev. B 48, 17006-17009 (1993). 
T. Arima and Y. Tokura, "Optical study of electronic 



structure in perovskite-type RMO3 {R = La, Y; M = 
Sc, Ti, V, Cr, Mn, Fe, Co', Ni, Cu)," J. Phys. Soc. Jpn. 64, 
2488-2501 (1995). 

I. S. Smirnova, "Normal modes of the LaMnOa Pnma 
phase: comparison with La2Cu04 Cmna phase," Physica 
B 262, 247-261 (1999). 

R. Kajimoto, H. Mochizuki, H. Yoshizawa, H. Shintani, T. 
Kimura, and Y. Tokura, "i?-dependence of spin exchange 
interactions in fJMnOa {R — rare-earth ions)," J. Phys. 
Soc. Jpn. 74, 2430-2433 (2005). 

T. Kimura, S. Ishihara, H. Shintani, T. Arima, K. T. Taka- 
hashi, K. Ishizaka, and Y. Tokura, "Distorted perovskite 
with configuration as a frusrated spin system," Phys. 
Rev. B 68, 060403(R) (2003). 

T. Goto, Y. Yamasaki, H. Watanabe, T. Kimura, and Y. 
Tokura, "Anticorrelation between ferromagnetism and fer- 
roelectricity in perovskite manganites," Phys. Rev. B 72, 
220403(R) (2005). 

"''^ J. Strempfer, B. Bohnenbuck, M. Mostovoy, N. Aliouane, 
D. N. Argyriou, F. Schrettle, J. Hemberger, A. Krimmel, 
and M. v. Zimmermann, "Absence of commensurate or- 
dering at the polarization flop transition in multiferroic 
DyMnOa," Phys. Rev. B 75, 212402 (2007). 
T. Arima, T. Goto, Y. Yamasaki, S. Miyasaka, K. Ishii, 
M. Tsubota, T. Inami, Y. Murakami, and Y. Tokura, 
"Magnetic-field-induced transition in the lattice modula- 
tion of colossal magnetoelectric GdMnOs and TbMnOs 
compounds," Phys. Rev. B 72, 100102(R) (2005). 

''^ A. A. Mukhin, V. Yu. Ivanov, V. D. Travkin, A. Pimenov, 

A. Loidl, and A. M. Balbashov, "Antiferromagnetic res- 
onance in the canted phase of Lai-^Sr^iMnOs: Experi- 
mental evidence against electronic phase separation," Eu- 
rophys. Lett. 49, 514-520 (2000). 

K. Hirota, N. Kaneko, A. Nishizawa, and Y. Endoh, "Two- 
dimensional planar ferromagnetic coupling in LaMnOg," J. 
Phys. Soc. Jpn. 65, 3736-3739 (1996). 
F. Kagawa, M. Mochizuki, Y. Onose, H. Murakawa, 
Y. Kaneko, N. Furukawa, and Y. Tokura, "Dynamics 
of multiferroic domain wall in spin-cycloidal ferroelectric 
DyMnOs," Phys. Rev. Lett. 102, 057604 (2009). 
''^ T. Moriya, "Far infrared absorption by two magnon excita- 
tions in antiferromagnets," J. Phys. Soc. Jpn. 21, 926-932 
(1966). 

B. R. Cooper, R. J. Elliott, S. J. Nettel, and H. Suhl, "The- 
ory of magnetic resonance in the heavy rare-earth metals," 
Phys. Rev 127, 57-68 (1962). 

®^ C. Jia, S. Onoda, N. Nagaosa, and J. H. Han, "Micro- 
scopic theory of spin-polarization coupling in multiferroics 
transition metal oxides," Phys. Rev. B 76, 144424 (2007). 



